Computational Analysis of the G-III Laminar Flow Glove 
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Under NASA’s Environmentally Responsible Aviation Project, flight experiments are 
planned with the primary objective of demonstrating the Discrete Roughness Elements 
(DRE) technology for passive laminar flow control at chord Reynolds numbers relevant to 
transport aircraft. In this paper, we present a preliminary computational assessment of the 
Gulfstream-III (G-III) aircraft wing-glove designed to attain natural laminar flow for the 
leading-edge sweep angle of 34.6°. Analysis for a flight Mach number of 0.75 shows that it 
should be possible to achieve natural laminar flow for twice the transition Reynolds number 
ever achieved at this sweep angle. However, the wing-glove needs to be redesigned to 
effectively demonstrate passive laminar flow control using DREs. As a by-product of the 
computational assessment, effect of surface curvature on stationary crossflow disturbances is 
found to be strongly stabilizing for the current design, and it is suggested that convex surface 
curvature could be used as a control parameter for natural laminar flow design, provided 
transition occurs via stationary crossflow disturbances. 

Nomenclature 

= wing/airfoil chord length 

= lift coefficient 

= pressure coefficient 

= disturbance phase speed 

= disturbance frequency 

= flight altitude 

= free stream Mach number 

= N-factor 

= Reynolds number based on wing chord 

= Reynolds number based on boundary layer momentum thickness 
= surface distance 

= adiabatic wall temperature 

= turbulence level 

= wall temperature 

= chordwise coordinate 

= spanwise coordinate 

= surface curvature 

= growth rate of disturbance 

= boundary layer momentum thickness 
= spanwise wavelength 

= wing sweep angle 

= disturbance wave angle 
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I. Introduction 


S Kin friction constitutes about 50 percent of the drag budget of a business jet or a long haul transport aircraft, 1 ' 3 
as depicted in Fig. 1 . Therefore, reduction in skin friction has the potential of yielding substantial savings in fuel 
burn. According to a system study, 3 10 percent reduction in skin-friction drag would result in about 9 percent fuel 
savings for a 4000 nm transport aircraft. Since laminar skin friction is much less than the turbulent value, the subject 
of maintaining laminar flow for high Reynolds number flows (i.e., delay of boundary-layer transition) has been of 
interest for well over half a century; see Refs. 4 ' 8 for reviews of the subject. 

In low disturbance environments, Tollmien-Schlichting (TS) instability is responsible for transition in two- 
dimensional (2D) boundary layers that can be controlled by favorable pressure gradients, wall suction and wall heat 
transfer (cooling in air; heating in water). In three-dimensional (3D) swept-wing boundary layers, crossflow 
instability dominates the transition process, which often manifests itself in the form of stationary corotating 
streamwise vortices that originate at minute roughness sites 9 ; and these vortices breakdown via a high frequency 
secondary instability mechanism. 10 ' 14 Crossflow instability can be controlled by using wall suction and there have 
been a number of experimental campaigns where laminar flow control (LFC) using suction has been demonstrated 
in subsonic 15 " 17 as well as in supersonic 18 flight. 

A three-dimensional boundary layer is subject to crossflow instability owing to an inflection point in the 
crossflow velocity component (i.e., the component in the direction normal to the inviscid streamlines). Minimizing 
this velocity component would weaken crossflow instability and yield extended runs of laminar flow. For a given 
leading-edge sweep angle and unit Reynolds number, the only way to reduce the crossflow velocity component is to 
minimize the pressure-gradient but that would bring TS instability into play. Thus, swept-wing Natural Laminar 
Flow (NLF) design is based upon striking a delicate balance by tailoring the pressure distribution such that the 
crossflow instability is reduced sufficiently while keeping TS growth at subcritical levels with respect to transition. 
One can only maintain this balance for a limited range of chord Reynolds numbers and LFC via suction must be 
employed if laminar flow is desired at high Reynolds numbers associated with large transport aircraft. 

The wavelength of crossflow disturbances scales with the boundary-layer thickness, with small wavelength 
disturbances growing first (i.e., at small distance from the wing leading-edge) followed by decay downstream. The 
longer wave length disturbances become unstable later but remain unstable over longer streamwise distances and 
eventually breakdown via secondary instabilities that lead to laminar-turbulent transition. Saric et al. 19 discovered 
that forcing small wavelength disturbances (i.e., the “control'’ mode) at relatively high amplitudes changes the 
boundary-layer mean flow such that the growth of more dominant larger wave length disturbances (i.e., the “target” 
mode) is pushed downstream, thus delaying transition. This concept has been studied by using nonlinear Parabolized 
Stability Equations (PSE) 20 and Direct Numerical Simulations (DNS) 21 and both these analyses support the 
soundness of the fundamental concept. In practice, the control mode is introduced via discrete spanwise periodic 
roughness elements placed near the neutral point of the instability diagram. These control disturbances only 
introduce harmonics and no subharmonics, and therefore do not directly feed into the larger wavelength target 
modes. This strategy for controlling crossflow disturbances is known as Discrete Roughness Elements (DRE) 
technology and it has the potential of increasing the range of applicability of NLF to higher chord Reynolds numbers. 

Saric et al. 22 " 23 reported a flight experiment using Cessna 0-2 aircraft to demonstrate the DRE concept. The test 
article consisted of a 30° swept blade that was mounted vertically under one of the wings of the aircraft. Maximum 
transition Reynolds number obtained on the test article was about 6.4 million when the leading edge was polished 
with surface finish quoted as 0.3 pm root mean square (rms) and 2.2 pm average peak-to-peak. Flere, transition 
Reynolds number was limited by pressure minimum (i.e., adverse pressure gradient) because of the relatively low 
chord Reynolds number of 8 million. In other words, increased extent of NLF could be had for a larger test article. 
When the leading edge was painted to yield a deteriorated surface finish (quoted as 1 pm rms, with 3.8 pm peak-to- 
peak), transition Reynolds number dropped to about 2.4 million owing to the adverse effect of surface roughness on 
transition; with the application of DRE, transition Reynolds number increased to about 4.8 million for the chord 
Reynolds number of 8 million. This is the highest Reynolds number at which DRE technology has been 
demonstrated to delay crossflow-induced transition. Clearly, natural laminar flow could be maintained at such 
Reynolds numbers without the use of DRE and there is need to demonstrate the technology at higher chord 
Reynolds numbers, Mach numbers and lift-coefficients of relevance to transport aircraft. 
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NASA's Environmentally Responsible Aviation (ERA) Project has selected DRE as one of the laminar flow 
control technologies for further evaluation and development for potential application to transport aircraft. A 
collaborative effort between Dryden Flight Research Center, Langley Research Center and Texas A&M University 
(TAMU) has been initiated whereby flight tests are being planned to increase the transition Reynolds number based 
on NLF as much as reasonably possible and then increase that Reynolds number even further with the use of DRE. 
An additional objective of the experiment is to relax the surface finish requirements for swept-wing NLF using DRE 
as suggested by the results of Ref 23. 

The aircraft selected for the flight experiments is Dryden’ s Gulfstream-III (G-III, hereinafter G-3) aircraft for 
which a wing glove has been designed by TAMU. 24 The glove leading-edge sweep is 34.6° with a maximum 
possible chord Reynolds number approaching 30 x 10 6 . The results in this paper provide a computational assessment 
of the glove using the unstructured-grid Navier-Stokes code, FUN3D 25 and stability analysis code, LASTRAC." 6 
Based on computational analyses performed at Langley and TAMU, the wing-glove is being redesigned, but the 
results reported herein use the glove geometry that is currently available. 

In this paper, we first describe the objectives of the DRE Laminar Flow Glove Experiment (DRE LFGE), also 
known as Subsonic Aircraft Roughness Glove Experiment (SARGE), in the context of previous NLF experiments. 
We then briefly describe the computational tools and present the results obtained to-date, while pointing out that the 
present analysis is only a snapshot of an ongoing effort as the glove geometry will change before it is manufactured 
and installed on the aircraft. The role of surface curvature in crossflow-induced transition is also studied. 

II. Research Background 


A. SARGE Objectives 

The goal set by the Environmentally Responsible Aviation Project is to demonstrate Discrete Roughness 
Elements technology for leading-edge sweep in excess of 30° and chord Reynolds numbers that are characteristics of 
mid-range transport aircraft, i.e.. Re c up to 30 x 10 6 . The Mach number range of interest is 0.66 to 0.75, with cruise 
section Cl of approximately 0.5. NASA Dryden's Gulfstream-III (G-3) aircraft was selected for the Discrete 
Roughness Elements Laminar Flow Glove Experiment (DRE LFGE), also known as the Subsonic Aircraft 
Roughness Glove Experiment (SARGE). The G-3 aircraft has a wing leading-edge sweep of 31.7°. The objectives of 
the planned flight tests are: 

• Demonstrate NLF for glove Re c > 15x10°, with laminar flow up to (x/c) tr = 0.6. Therefore, the minimum 
required transition Reynolds number on the suction side of the glove is 9 million for the prescribed sweep 
angle. 

• Demonstrate effectiveness of DRE in delaying transition by at least 50% at Re c > 22x1 0 6 . This means that if 
natural transition for a given surface finish occurs at (x/c) tr = 0.4, DRE will delay it to (x/c) tr = 0.6. 

• Perform experiments with polished and painted leading-edges to demonstrate the effectiveness of DRE in 
relaxing surface quality requirements for NLF at high Reynolds numbers. 

Two rows of pressure ports will be provided to measure surface pressure and thermocouples will be used to 
measure wall temperature. Transition location on the suction side will be inferred by using an Infra-Red camera and 
hot films. Laminar flow extent will be measured on the pressure side by using thermocouples, but maximizing 
laminar flow on this side is not one of the primary objectives of the experiment. The experiment should also yield 
valuable data for the validation of computational tools. 

B. Previous NLF Experiments: 

Several Natural Laminar Flow (NLF) experiments have been conducted over the years in wind tunnels and 
during flight experiments; and some of the results are presented in Fig. 2 where maximum transition Reynolds 
numbers are plotted as a function of the sweep angle. Here, transition Reynolds number is based on free stream 
conditions, rather than local conditions. Most of the data points are taken from Ref. 27 and plotted here for 
comparison with the intended range of the planned G-3 flight experiment. These include results from the B-18 flight 
test, 28 King Kobra flight test, 29 the Boeing 757 NLF glove, 30 F-l 11/TACT NLF glove flight experiment, 31 F-14 
variable sweep transition flight experiment, 32 low-speed experiments conducted in the Ames 12-Foot Tunnel, 33 and 
Langley Transonic Pressure Tunnel (TPT). 34 More recent sets of data include the White-Knight experiment 35 


3 

American Institute of Aeronautics and Astronautics 



conducted by Northrop-Grumman and the Cessna 0-2 flight experiment referred to above. Both these experiments 
were conducted on 30° swept blades installed under the aircraft. 

Largest values of transition Reynolds number were obtained for the F-14 flight experiment, where “... a 
maximum transition Reynolds numbers of about 17.6 million was obtained at 15 degrees of sweep, 13.5 million at 
20 degrees, and 12 million at 25 degrees. Beyond 25 degrees of sweep, maximum transition Reynolds number 
decreased rapidly to about 5 million at 35 degrees of sweep”. 27 It was noted that maximum transition Reynolds 
numbers were limited by either adverse pressure gradients (i.e., TS dominated transition) or shock wave location. 
Thus, for a larger wing chord Reynolds number, it should be possible to obtain even higher values of crossflow- 
induced transition Reynolds numbers for an NLF wing. . 

Transition Reynolds numbers for flight experiments conducted in Europe (e.g., ATT AS 36 and F-100 Fokker 37 ) 
are not available in the open literature. Flowever, Ref. 37 notes that in the latter experimental campaign it was 
possible to achieve laminar flow for Reynolds numbers up to 25 million, for leading edge sweep angles below 23°. It 
is reasonable to assume that the highest value of transition Reynolds number of 25 x 10 6 was associated with the 
lowest sweep angle of 17° for the F-100 experiment. This is to be compared with the corresponding value of about 
16 x 10 6 for the F-14 experiment. If the above inference from Ref. 37 is correct, the highest reported transition 
Reynolds number achieved in an NLF flight experiment is from the F-100 Fokker experiment. Flowever, the authors 
have not been able to confirm this with the principals involved in the experiments. 

Transition Reynolds number decreases with increasing sweep angle due to the destabilizing effect of crossflow. 
Data points from two recently conducted flight experiments (White-Knight 35 and Cessna 0-2 23 ) have also been 
included in Fig. 2. The test article in the White-Knight experiment also included a 30° swept blade attached to the 
aircraft. Maximum transition Reynolds number obtained in this case was about 5.6 x 10 6 and was limited by the 
pressure minimum (i.e., adverse pressure gradient) because of the relatively low chord Reynolds number of 7. 1 x 10 6 . 
Results from the Cessna 0-2 have been discussed before. The expected range of transition Reynolds number of the 
SARGE is also shown in Fig. 2, where the lower bound comes from the NLF requirement and the upper bound is 
derived from the maximum Reynolds number attainable for the G-3 glove under the assumption that laminar flow 
can be achieved up to 60% chord using DRE. The minimum expected Reynolds number of 9 x 10 6 is 80% higher 
than that achieved in the F-14 flight experiment at this sweep angle of about 35 degrees. 

C. Glove Design 

Glove design was carried out by Texas A&M University and is discussed in Ref. 24. The glove as installed on 
the G-3 aircraft wing is shown in Fig. 3(a-b). It has a leading-edge sweep of 34.6°, is 6 ft wide and has a streamwise 
chord length of 14.5 ft. The glove is smoothly faired into the wing on inboard, outboard and downstream sides. With 
wing-glove thickness ( t/c ) and section-lift coefficient (Cl) as constraints, design proceeded to place the pressure 
minimum as far aft as possible so that the boundary-layer is subcritical to TS-induced transition. Linear stability 
results were used to guide the design for crossflow control. The leading-edge radius was limited by the requirement 
that Reg < 100 to assure that the attachment-line boundary layer is subcritical to any span wise contamination, 
although this requirement could have been relaxed owing to the inboard leading-edge fairing that would act as a 
Gaster bump. The requirement for clearance between the glove and the wing resulted in a thicker and longer glove 
(see Fig. 3(c)) than it might have been otherwise. Flowever, this inadvertently made the design more interesting from 
the view point of the effect of surface curvature, which we will address towards the end of the paper. 

III. Computational Tools 

A. FUN3D 

FUN3D 25 ' 38 is a finite-volume Navier-Stokes solver in which the flow variables are stored at the vertices or nodes 
of the mesh. FUN3D solves the equations on mixed element grids, including tetrahedra, pyramids, prisms, and 
hexahedra and also has a two-dimensional path for triangular/quadrilateral grids. At interfaces delimiting 
neighboring control volumes, the inviscid fluxes are computed using an approximate Riemann solver based on the 
values on either side of the interface. Several convective flux schemes are available in FUN3D. The most common 
scheme for subsonic and transonic flows is Roe’s flux difference splitting, 39 which is used in the current study. For 
second-order accuracy, interface values are obtained by extrapolation of the control volume centroidal values, based 
on gradients computed at the mesh vertices using an unweighted least-squares technique. 

For tetrahedral meshes, the full viscous fluxes are discretized using a finite-volume formulation in which the 
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required velocity gradients on the dual faces are computed using the Green-Gauss theorem. On tetrahedral meshes 
this is equivalent to a Galerkin type approximation. For non-tetrahedral meshes, the same Green-Gauss approach can 
lead to odd-even decoupling. A pure edge-based approach can be used to circumvent the odd-even decoupling issue, 
but yields only approximate viscous terms. Thus for non-tetrahedral meshes, the edge-based gradients are combined 
with Green-Gauss gradients, which improves the h-ellipticity of the operator, and allows the complete viscous 
stresses to be evaluated. 40 This formulation results in a discretization of the full Navier-Stokes equations. 

The solution at each time-step is updated with a backwards Euler time-integration scheme. At each time step, the 
linear system of equations is approximately solved with either a multi-color point-implicit procedure or an implicit- 
line relaxation scheme. 38 Local time-step scaling is employed to accelerate convergence to steady state. For 
turbulent flows, several models are available within FUN3D, including the loosely coupled Spalart-Allmaras (S-A) 
model. 41 

B. Unstructured Grid 

The unstructured-grid tool VGRID 42 was employed to generate a tetrahedral grid for the G-3 configuration with 
the glove. Only one-half of the aircraft was gridded as flow symmetry was assumed along the centerline of the 
fuselage. Since the region of particular interest for this study is the mid-span region of the wing, the tail was not 
simulated and some minor details of the engine nacelle geometry were ignored. Powered engine conditions as 
provided by Ref. 43 were used for computation of the aircraft flowfield. A mixed-element grid consisting of 
tetrahedrons and prisms was used in the viscous flow computations, where the prisms were generated by 
recombining the tetrahedral grid in the advancing layers. 

C. LASTRAC 

Stability calculations have been performed by the Langley Stability and Transition Analysis Code, LASTRAC." 6 
The LASTRAC software is a suite of general purpose, stability and transition analysis tools that encompass a broad 
range of the state-of-the-art numerical methods based on linear stability theory (LST), linear and nonlinear 
parabolized stability equations (PSE), and adjoint PSE. Stability calculations and N-factor correlations can be 
performed for both 2D/axisymmetric and general 3D configurations. 44 Nonlinear PSE modules in LASTRAC have 
been used to study the effect of discrete roughness elements on crossflow transition. 45 The governing equations in 
LASTRAC account for surface curvature and nonparallel mean flow. 

IV. Results and Discussion 

A. Pressure Coefficients 

Calculations were first performed for the case with no glove to compare with available experimental data to 
validate the CFD results. Good agreement (not shown) was achieved for the pressure distribution and lift coefficient. 
Computations were then carried out for the G-3 aircraft with the gloved swept-wing. The glove has a slightly larger 
leading-edge sweep of about 34.6° compared to 31.7° for the rest of the wing. The flow conditions for the present 
computations are: M = 0.75, Re c = 22 x 10 6 . H = 38,840 ft, and angle of attack ( AoA) = 3.5°. The unstructured-grid 
generated for this case is depicted in Fig. 4. The grid has a total of about 27 million nodes. The solution is obtained 
in two ways. First, turbulent flow is computed everywhere using the S-A model and the solution in this case 
converges rapidly. Second, turbulent flow was assumed everywhere except on the glove up to just upstream of the 
shock and this solution also converged. Results for the two computations are compared in Fig. 5. Pressure 
coefficients are essentially the same in the two cases except near the shock where the partly laminar computation 
shows a kink due to the abrupt turn-on of the turbulence model. We note that the boundary-layer grid used here is 
not sufficiently fine for this Navier-Stokes solution to be used directly for the stability analysis. Therefore, such 
analysis will be based on boundary-layer computations using the Cp distributions derived from the Navier-Stokes 
solution and for that purpose either a turbulent or laminar CFD computation for the glove could be used. 

Figure 6 gives pressure distributions for three selected span stations (y = 204 in, 234 in and 264 in) on the glove. 
Only the suction side results are of particular interest here. The Cp results generally exhibit favorable gradients as 
desired. Along the mid-span station on the glove, y = 234”, the pressure gradient is favorable all the way to the 
shock and should result in crossflow dominated transition. Along the inboard glove station y = 204”, there is a strong 
favorable pressure gradient up to x/c = 0.14 which would result in stronger crossflow instability at this inboard 
location as compared to the y = 234” station. Flowever, the nearly flat roof top distribution between x/c = 0. 14 and 
0.4 is expected to bring TS instability. At the outboard station y = 264”, the Cp distribution is nearly flat between 
x/c = .02 and .2; therefore, TS instability should dominate in this region, followed by stronger crossflow growth 
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downstream of x/c = .2 where stronger favorable gradients are present. These observations will be confirmed by 
performing linear stability computations next. 

B. Linear Stability Analysis of the G-3 Glove Boundary-Layer 

Streamwise pressure distributions at the three spanwise stations on the glove were extracted from the CFD 
solution and used for computing mean flow using the boundary layer code BLSTA. 46 The boundary layer profiles 
thus computed were analyzed using the LASTRAC code under the quasi-parallel flow assumption. This mode of 
analysis is commonly referred to as Linear Stability Theory (LST). Three-dimensional boundary layers are unstable 
to both stationary and traveling crossflow disturbances, with traveling disturbances having much higher growth rates 
based on linear stability analysis. The question of whether stationary or traveling disturbances are relevant in the 
transition process primarily depends upon the receptivity of these disturbances to external forcing. Surface 
roughness induces stationary disturbances, while free-stream turbulence induces traveling disturbances. The effect 
of free stream turbulence on crossflow-induced transition was studied by Deyhle and Bippes 47 in wind tunnels, 
where it was found that stationary disturbances dominate transition at low turbulence levels, while traveling 
instability waves become dominant when turbulence levels exceed 0.2 percent ( Tu > 0.2). Saric et al. 23 have argued 
that since the turbulence level in flight is low, only stationary disturbances are relevant in flight transition making 
the boundary-layer subject to DRE control. Therefore, the present analysis is limited to stationary disturbances. It is 
important to note that in cases where traveling disturbances dominate the transition process, DRE might not be 
relevant in delaying transition. 

For stationary (f= 0 Hz) crossflow instability analysis, computations were performed for disturbances of various 
spanwise wavelengths which were held fixed as the N factor computations proceeded in the downstream direction. 
Here, N is given as 


N = / a ds (1) 

where a is the disturbance growth rate computed using stability analysis and ,v represents the surface distance along 
the path of integration, which is either in the chordwise direction or along the inviscid streamlines. Disturbances 
grow by a factor of e N prior to transition, where the value of N is empirically determined. 

The results for disturbance wavelengths of 3 mm to 18 mm are shown in Fig 7, which confirm our conjecture 
based on the Cp distributions. At the inboard location, crossflow N factor increases rapidly and a value of 9 is 
reached at x/c = 0.14. At the mid-span of the glove, N factor growth is more gradual and an N of 9 is reached at 
about x/c = .44. For the outboard location, N factors are much smaller initially but stronger growth downstream of 
x/c = .25 yields N = 9 at x/c = .44 as for the mid-span location. We use the value N = 9 for illustrative purposes and 
do not imply that this would necessarily correlate with transition in this case. The value of N at transition depends 
upon surface finish near the leading edge. Saric et al. 23 noted that N reached a value of 14 for the polished leading 
edge in the Cessna 0-2 flight experiment. For the G-3 glove, this value is reached at x/c = .52, .58 and .56 for the 
inboard, mid-span and outboard locations, respectively. 

Computations are next performed for TS disturbances of various frequencies ranging from 1,500 to 6,000 Hz; 
the results are presented in Fig. 8. These results are obtained by holding wave angle fixed with respect to the 
inviscid streamline. The value of wave angle is selected to be 0° in this case, i.e., the wave vector is assumed to align 
with the streamline. The results show that the N factor for TS waves is the highest (N = 6) at the inboard location 
followed by the outboard location (N = 5.5) with mid-span location having the lowest value (N = 4). These values 
are again consistent with our expectations from the Cp distributions. The effect of TS wave orientation is studied in 
Fig. 9 where we show results for wave angles of -20°, 0°, 20° and 30° for the mid-span location. The N factor for TS 
waves increases with wave angle, reaching 5.5 at 30° as compared to 4 for the 0° wave angle case. Some comments 
on the TS N factor computation strategies will be made in the next section. 

The above computations have been made assuming adiabatic wall conditions. During a short duration test flight, 
the glove temperature would start out from being at atmospheric value at sea level. As the aircraft approaches the 
design altitude, the surrounding temperature drops much more rapidly than the glove temperature. Therefore, before 
an equilibrium can be established, conditions of T w > T adb can exist. It is well known that TS instability is very 
sensitive to wall heat transfer and, in particular, experimental data in Ref. 48 shows that transition Reynolds number 
varies as (T w /Tadb)" 7 . In order to assess the wall temperature effect for the glove boundary layer, calculations were 
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made by assuming T w = l.lTadb for the mid-span location assuming wave orientation of 30°. Results in Fig. 10 show 
that the N factor reaches a value of 8 for this temperature condition as compared to 5.5 for the adiabatic case in Fig. 
9. It is important that the NLF experiment be performed only after the glove temperature has dropped sufficiently 
close to the adiabatic wall temperature. This can be accomplished via a ‘cold soak' where the aircraft flies at a 
higher altitude (i.e., colder atmospheric temperature) before descending to the prescribed test altitude. 

Based on the significant spanwise variations of Cp on the G-3 glove noted above, and the associated stability 
results, it is recommended that the glove be redesigned to minimize the spanwise variations. This redesign is 
currently underway. 

C. PSE Analysis 

As noted above, swept-wing transition involves nonlinear growth of crossflow disturbances followed by 
breakdown via high frequency secondary instabilities 11 and this evolution is readily computed by nonlinear 
Parabolized Stability Equations (PSE) and secondary instability analyses. 45 A crucial piece of missing information is 
the initial disturbance level induced by roughness. Carpenter et al. 49 have studied receptivity of incompressible 
swept-wing boundary layer to surface roughness using Linearized Navier-Stokes (LNS) equations. For the Cessna 
0-2 case, they were able to demonstrate that disturbance amplitude internalized in the boundary-layer differed by 
two orders of magnitude (approximately e 5 ) between surface roughness distributions modeled after surface height 
measurements for painted and polished surfaces, which is roughly consistent with the observation of Saric et al.~ 3 
that transition moved forward to about 30% chord when the leading-edge was painted as compared to the transition 
location at 80% chord for the polished leading-edge. We note that Ref. 49, after detrending the raw data of Ref. 23, 
the computed rms surface roughness height was 0.09 pm and 6.5 pm for the polished and painted surfaces, 
respectively, as compared to the raw rms values of .3 pm and 1 pm quoted in Ref. 23. 

Receptivity studies require spectral information about surface finish, which will not be available for the G-3 
glove until it is designed, manufactured and measured. In the absence of such knowledge, the only alternative is to 
assume an initial disturbance level and then parametrically study its impact on secondary instability breakdown. 
Such an analysis was conducted by Ref. 45 for an infinite swept-wing to investigate DRE effectiveness. Their 
analysis was performed for a swept airfoil glove 50 intended for a Gulfstream-II (G-2) aircraft wing, a hypothetical 
predecessor to the current G-3 effort. Computations were performed for a Mach number of 0.75, chord Reynolds 
number of 17 x 10 6 and wing leading-edge sweep of 30°. Nonlinear PSE analysis was used to first compute the 
development of stationary crossflow vortices and then secondary instability analysis was performed with a view to 
assess the feasibility of roughness based control methodology at high Reynolds numbers. It was noted that a 
broad spectrum of stationary crossflow disturbances amplify and. while it may be possible to control a specific 
target mode using Discrete Roughness Elements (DREs), nonlinear interaction between control and target modes 
may yield strong amplification of the difference mode that could have an adverse impact on the transition delay 
using spanwise periodic roughness elements.” Since these conclusions were drawn based on assumed initial 
amplitudes of disturbances, they can be characterized as qualitative at best. Definitive conclusions can only be 
drawn when initial modal amplitudes are available from a careful receptivity analysis for a given surface finish. 

A Compressible LNS code is under-development but is currently not available to carry out receptivity 
computations for the G-3, similar to those in Carpenter et al. 49 In the meantime, it will be instructive to carry out 
PSE analysis as a risk reduction exercise, particularly since the planned range of flight Reynolds numbers is even 
higher than the G-2 case of Ref. 45. Nonlinear PSE analysis using the LASTRAC code can be performed only for 
an infinite swept-wing and therefore rest of the analysis will be performed under this assumption. The G-3 glove 
geometry (Fig. 3) was developed from an infinite-swept airfoil design and we compare the Cp distribution for that 
airfoil (i.e., target Cp for the glove design) with the mid-span Cp for the 3D glove in Fig. 11. The two distributions 
are similar near the leading-edge but there is a large change in slope of the glove Cp ahead of the shock, which gives 
rise to larger N factors in that region noted earlier in Fig 7. N factor computations based on linear stability analysis 
for the infinite-swept airfoil are shown in Fig. 12 where maximum N factors for the 22 x 10 6 case reach about 10.6 
at x/c = 0.6 as compared to 14 in Fig. 7 (b). According to Saric’s LST-based e 14 criterion 23 for polished leading- 
edges, transition should take place on the glove at x/c = 0.58 for the current Reynolds number of 22 x 10 6 , but the 
boundary-layer would remain laminar if a redesigned glove can maintain the target infinite-swept wing pressure 
distribution shown in Fig. 11. Figure 12 also contains N factor results for Re c = 17 x 10 6 and 27.5 x 10 6 and 
transition will not occur in the highest Reynolds number case for the polished leading-edge based on the e 14 LST 
criterion. Transition for the target Cp distribution should occur at higher chord Reynolds numbers or with painted 
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(roughened) leading-edges. In order to achieve transition at Re c = 22 x 10 6 with polished leading-edge, the slope of 
the Cp curve must be increased and this should also help stabilize the TS waves. 

The first step in nonlinear analysis is to perform linear PSE computations to determine control and target modes, 
as was done in Ref. 45 Their results presented in Fig. 13, for example, show that stationary modes over a broad 
range of span wise wavelengths (4 mm < X <17 mm) attain a peak N factor greater than 8, based on linear PSE 
analysis of the G-2 airfoil. The 8 mm wavelength achieved the highest N factor of about 12.3, and this wavelength 
was selected as the target mode with the 2.67 nun wavelength acting as the control mode in their nonlinear analysis. 
It would be instructive to perform linear PSE analysis for the G-3 design airfoil to determine the target and control 
modes for DRE application. These results are presented in Fig. 14 and it is clear that the highest N factor in this case 
is substantially lower than the LST result shown in Fig. 12. It can thus be argued that if stationary modes constitute 
the dominant route to breakdown, transition will not take place for this case with Re c = 22 x 10 6 as the highest N 
factor is only about 4.6. The PSE N factor is lower by almost 6 as compared to the results based on LST, which is an 
unexpected result and it is important to investigate this further. Therefore, rest of the paper will focus on this 
dramatic drop in N factor and a nonlinear PSE analysis is deferred to a future date. 

Parabolized Stability Equations account for nonparallel and surface curvature effects, while these effects are 
ignored in LST. References 51 " 2 showed that nonparallel effect is usually destabilizing and surface curvature is 
stabilizing as compared to computations based on the quasi-parallel assumption of LST. Comparison of the results in 
Fig. 12(b) and Fig. 14 implies that the surface curvature effect is much stronger than the nonparallel effect in this 
case. We study this in Fig. 15 where we present the results for X = 9mm using (a) quasi -parallel theory with no 
curvature, (b) quasi-parallel theory with curvature, (c) PSE with no curvature (i.e., nonparallel effect only) and (d) 
PSE with curvature. The last option is the only consistent approach in that it accounts for all the relevant effects of 
the same order. The N factor using PSE drops from 14.5 in the no-curvature case to only about 4.1 when curvature 
terms are included. Figure 16 presents the results for the G-2 airfoil with and without curvature using X = 8mm and 
in this case the N factor drops from about 17 (PSE, no curvature) to 12.3 (PSE, with curvature). While the two (G-2 
and G-3) Cp distributions are different, it is still remarkable that the maximum N factor (for X = 12mm) is only 4.6 
for the G-3 airfoil at Re c = 22 x 10 6 as compared to 12.3 for the G-2 airfoil at the lower Re c = 17 x 10 6 . Thus, the 
drop in maximum N factor due to curvature is 4.7 and 9.9 for the G-2 and G-3 airfoils, respectively; the stabilizing 
effect is much stronger in the latter case. Since the maximum N factor difference between LST and PSE is about 1 .7 
and 6 for the G-2 and G-3 airfoils, respectively, use of the LST approach in design computations is justified in the 
former case but will introduce large errors in the latter case. 

The centrifugal term in the normal momentum equation that induces Gortler instability 53 in boundary-layers over 
concave surfaces is also the term that induces the stabilizing effect of convex curvature on crossflow vortices. 51 
Therefore, Gortler number should provide a measure of the effect of convex curvature on crossflow instability as it 
does for centrifugal instability on concave surfaces. In order to understand the disparate effect of curvature for the 
G-2 and G-3 airfoils, we plot airfoil geometries and Gortler number for the two cases in Fig. 17. Here, Gortler 
number is defined as 


G g = Re g yf6ic (2) 

where 0 is the chordwise boundary-layer momentum thickness and k is surface curvature. The comparison in Fig. 17 
clearly shows that the Gortler number is much higher for the G-3 airfoil; and this explains the stronger stabilizing 
effect in this case. Figure 17 also gives Gortler number distributions for the G-3 case with Re c = 17 x 10 6 and 27.5 x 
10 6 . The oscillations in the Gortler number for the G-3 are the direct result of oscillations in surface curvature and 
this induces oscillations in the growth rate noticed in Fig. 15. 

N factor results for the three Reynolds numbers (Re c = 17 x 10 6 , 22 x 10 6 and 27.5 x 10 6 ) are compared in Fig. 18. 
The peak N factors for the three cases are 3.3, 4.6 and 6.6, respectively. The value of 3.3 is to be compared with N = 
12.3 for the G-2 case at the same Reynolds number of 17 x 10 6 . One thing is clear from this analysis: for the lowest 
Reynolds number, the boundary-layer will remain laminar up to x/c = 0.6 yielding a transition Reynolds number of 
10.2 x 10 6 . This is almost twice the Reynolds number achieved in the F-14 NLF flight experiment at the same sweep 
angle. 
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If transition is induced by stationary crossflow vortices (a necessary condition for the applicability of DRE 
technology), then the above results suggest that transition will not occur for the two lowest Reynolds numbers and, 
therefore, the DRE control experiments are not relevant even for the Re c = 22 x 10 6 case. Crossflow-induced 
transition is likely for the highest Reynolds number case but perhaps at relatively large x/c with little scope for DRE 
control, except with roughened glove leading-edges. However, at this Reynolds number, TS N factors reach values 
in excess of 8 for adiabatic wall conditions and TS instability may become the dominant cause of transition. 
Therefore, it is important that the glove design be changed to allow steeper favorable pressure gradients to suppress 
TS waves. This will also result in higher stationary crossflow disturbance growth using PSE, which can then be 
suppressed or delayed via DRE in order to address the primary objective of the flight experiment. 

It is important to reiterate that the success of DRE rests upon the premise that traveling crossflow disturbances 
are not generated in the 3D boundary-layer because of the low turbulence levels in flight 23 . Based on the 
experimental data of Refs. 47 ' 54 , it was argued in Ref. 55 that transition is caused by stationary disturbances when Tu 
~ .02% and dominated by traveling instability waves when Tu ~ .15%. Masad and Abdelnaser 56 analyzed the data 
from the swept cylinder experiment of Ref. 57 and argued that transition was induced by traveling instability waves 
in this case where Tu ~ .15%. Their N factor analysis, which included nonparallel and surface curvature effects, 
yielded mean N values of 14.9 and 7.4 for traveling and stationary disturbances, respectively. Takagi and Itoh 5S 
noted that, as they continued to clean the surface of their swept cylinder, stationary disturbances gave way to 
traveling instability waves in their low disturbance wind tunnel experiment. It should be noted, however, that the 
swept cylinder used in Ref. 58 had an order of magnitude higher leading-edge radius than the G-3 airfoil and the 
associated thicker boundary-layer in that case would be much less sensitive to surface roughness. For the 757 HLFC 
flight experiment, transition location could be well-correlated with a traveling disturbance N factor of 9.6 calculated 
using PSE. This result does not necessarily imply that transition in the 757 experiment was caused by traveling 
instability waves as leading-edge suction distribution is known to induce a roughness effect, which would yield 
higher initial amplitude for stationary disturbances requiring smaller amplification by the boundary layer via linear 
mechanism. 

Traveling disturbances, if excited at a sufficiently high amplitude, grow much more than the stationary ones and 
in that case would dominate transition. Linear PSE results for traveling disturbances for Re c = 22 x 10 6 are presented 
in Fig. 19 using LST and PSE. First, we note that PSE N factor values are marginally lower than the LST value in 
this case. Second, the highest PSE N factor is about 15.6 as compared to only 4.6 in Fig. 14(b) for the stationary 
disturbances. Thus, PSE N factors (that include nonparallel and surface curvature effects) differ by almost 11 in this 
case for traveling and stationary disturbances. The difference drops to about 6 when LST is used in the analysis. It is 
not known if traveling disturbances could begin to dominate crossflow transition when N factors differ by about 1 1 . 
However, whether the difference in traveling and stationary N factors at small x/c is sufficiently large to overcome 
their disparate initial amplitudes is the only relevant question. This is because when disturbances grow to nonlinear 
amplitudes, results of linear analysis are no longer applicable and an N factor difference of 1 1 near the end of the 
domain of interest may be of no significance. 

As noted above, DRE design requires that the boundary-layer is subcritical to TS-induced transition. In 
computing TS N factors (Figs. 8-10), the engineering approach of fixing wave angle with respect to in viscid 
streamlines (cf.. Ref. 59) was employed. However, there is no physical basis for this approach and the only rational 
way to perform N-factor analysis for infinite span wings is to assume constant spanwise wavelengths as was done in 
the crossflow analysis above (see also Ref. 60). For the G-3 airfoil, TS wave computations have been carried out 
using LST for various disturbance frequencies and the frequency of 2000 Hz gave the highest N factor. Results of 
these computations for the most amplified frequency are shown in Fig. 20(a) for several wave angles ranging from 
0° to 50°. N factor continues to increase with wave angle as shown in the figure. The PSE results, with holding 
spanwise wavelengths fixed, are shown in Fig. 20(b) for the 2000 Hz frequency and several different disturbance 
wavelengths. Stability characteristics (N factor, disturbance growth rate, wave angle and phase speed) derived from 
these PSE computations are presented in Fig. 21 for two locations (x/c = .25 and .45). Wave angle increases with 
decreasing wave length, with highest N factor for a wave angle of about 70°. Instability for the selected frequency 
belongs to both the TS and traveling crossflow type. The disturbances with larger wavelengths belong to the TS 
family (viscous instability) while the disturbances with relatively smaller wavelengths are of traveling crossflow 
(in viscid instability) type. It is clear from the results that the switchover from one instability type to the other is not 
abrupt. For wave lengths X > 60 mm, phase speed is almost constant at about .4 when scaled with the free stream 
velocity, which is typical of TS waves. At small wave lengths, phase speed drops to a value of about .05, which is 
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typical of traveling crossflow disturbances. At intermediate wave lengths, the instability is of mixed type. Traveling 
disturbances are sensitive to free stream turbulence level and, at low levels of turbulence, TS waves are more 
sensitive to the acoustic disturbance field. Therefore, whether or not traveling crossflow and TS waves are relevant 
in the transition process depends upon free stream turbulence level at flight altitudes and the aircraft acoustic field. 

Definitive conclusions about transition mechanism in a 3D boundary-layer can only be drawn when information 
about external forcing as well as computational tools for receptivity prediction are available. It should now be 
possible to perform such an analysis to determine the effect of surface roughness (see Refs. 45 ' 49 ). However, we do 
not have sufficient information about turbulence characteristics (turbulence level and spectra) at flight altitudes and, 
therefore, quantitative prediction of receptivity of the G-3 glove boundary-layer to free-stream disturbance 
environment is not currently feasible, without which definitive conclusions cannot be drawn to completely rule out 
the role of traveling disturbances in flight and hence the need for physical experiments. If DRE control is 
demonstrated in flight at the selected Reynolds numbers, then that will provide an indirect evidence that traveling 
disturbances are not relevant in flight transition for transport aircraft even though they amplify much more within 
the framework of linear PSE. If the DRE control fails, then either the premise that only stationary crossflow 
disturbances are relevant in flight transition at all Reynolds numbers would be false or that the boundary-layer was 
not subcritical to TS-induced transition. 


V. Conclusions 

Computational analysis has been performed for the laminar flow glove designed for the Gulftream-III flight 
experiment planed under NASA's Environmentally Responsible Aviation Project to evaluate the effectiveness of the 
Discrete Roughness Elements (DRE) technology at high Reynolds numbers of relevance to transport aircraft. The 
Mach number, lift coefficient and sweep angle used in the analysis are 0.75, 0.5 and 34.6°, respectively. This 
analysis is limited to the suction side, which is the designated test surface for demonstrating passive laminar flow 
control using DREs. Computational tools included FUN3D~ 5 unstructured-grid Navier-Stokes code and LASTRAC 26 
for analyses based on Linear Stability Theory (LST) and Parabolized Stability Equations (PSE). The following 
conclusions can be drawn from the analysis: 

• For the current design, natural laminar flow should be achievable up to x/c = 0.6 at Re c = 17 x 10 6 yielding 
a transition Reynolds number of 10.2 x 10 6 , which is twice the Reynolds number ever achieved in an NLF 
experiment at the selected sweep angle of 34.6°. 

• There are significant spanwise variations in Cp; therefore, the glove needs to be redesigned to minimize 
these variations. 

• To demonstrate DRE technology, the new design needs to have stronger favorable pressure gradient in 
order to weaken TS growth, and to support stronger stationary crossflow instability that could then be 
controlled via DRE. 

• During flight testing, attention needs to be given to maintain adiabatic wall temperature for TS control. 

• As an unexpected outcome of the computational analysis, surface curvature was found to exert a very 
strong stabilizing influence on stationary crossflow disturbances for the present design. PSE N factors 
(including nonparallel and surface curvature effects) for stationary disturbances are smaller by about 6 
when compared to LST (quasi-parallel, without curvature) results for the three chord Reynolds numbers 
(Re c = 17 x 10 6 , 22 x 10 6 and 27.5 x 10 6 ) computed. For the intermediate Reynolds number case, the net N 
factor drop of 6 consisted of an increase in N factor of about 4 owing to the nonparallel effect and a 
decrease in N factor of about 10 owing to the surface curvature effect. 

• If in-flight transition indeed happens via stationary crossflow disturbances, then transition may also be 
delayed by using convex surface curvature as a control parameter. 

• In contrast with the above results for stationary disturbances, the PSE N factor for traveling crossflow 
disturbances was lower by only about 1 when compared to the values obtained using LST. The disparate 
effect of using PSE helps increase the difference in overall growth of traveling and stationary disturbances. 
PSE N factors for traveling disturbances are higher by almost 1 1 as compared to the values computed for 
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stationary disturbances; the difference drops to about 6 when LST is used (i.e., nonparallel and surface 
curvature effects are ignored) in the analysis. 

• The role, if any, of traveling crossflow disturbances during transition over the present glove design cannot 
be currently ascertained computationally due to the lack of knowledge (level and spectra) of free stream 
turbulence at the flight altitudes. 
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Figure 2. Maximum transition Reynolds numbers as a function of wing sweep for various NLF experiments. 
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(c) 

Figure 3. Gulfstream-III (G-3) aircraft and the glove, (a) Glove as installed on the G-3 wing, (b) Planform 
showing glove dimensions and fairings, (c) Cross-section of wing-glove. 




(a) (b) 

Figure 4. Tetrahederal surfacegrid used in the analysis, (a) Grid for the G-3 aircraft, (b) Grid in the glove region. 
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Figure 5. Computed Cp at the mid-span location (y = 234”) assuming turbulent flow everywhere and turbulent flow 

with laminar on the glove. 



Figure 6. Comparison of Cp distributions at three span wise stations on the glove, y =204”, 234”, and 264”. 
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Ca) y =204” 



(c) y =264” 

Figure 7. Stationary crossflow N factors computed using LST for several disturbance wavelengths. 
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Figure 8. LST N-factor results for TS waves computed for fixed wave angle (\|/) of 0° with respect to inviscid 

streamlines. 
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Figure 9. Effect of wave angles on TS N factors: (a) \) / = 0°, (b) \) / = 20°, (c) \| / = 30°, (d) \)/ = -20°. 
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Figure 10. TS N factors with T w = 1.1 T adb and \|/ = 30° 



(compare with Fig. 9 (c)). 


3D glove with that of the infinite-swept airfoil (target Cp). 
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Ca) Re c =17xl0 6 



(b) Re c =22xl0 6 



Figure 12. LST crossflow N factors for the infinite-swept airfoil. 
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Figure 13. PSE N factor results for the G-2 airfoil, Re c = 17 x 10 6 (from Ref. 45) 
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Figure 14. PSE N factor results for the G-3 airfoil, Re c = 22 x 10 6 (LST result for the most amplified wavelength 

shown for comparison). 




Figure 15. Nonparallel and surface curvature effects on disturbance growth rate and N factor with X = 9mm for G-3 

airfoil, Re c = 22 x 10 6 . (a) Growth rate, (b) N factor. 
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Figure 16. Nonparallel and surface curvature effects on N factors for the G-2 airfoil, Re c = 17 x 10 6 . 



Figure 17. Comparison of G-2 and G-3 airfoils: (a) chordwise geometry, (b) Gortler number. 
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Figure 18. PSE N-factor results for the three Reynolds numbers (for reference, LST results for most amplified wave 

length are also shown). 
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(a) (b) 

Figure 19. N factors for traveling instability waves (f= 1,000 Hz) of different spanwise wavelengths, (a) LST, (b) 

PSE. 
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Figure 20. N factors for a disturbance frequency of 2,000FIZ. (a) LST for different wave angles, (b) PSE for 

different spanwise wavelengths. 



Figure 21. PSE growth rate (a), N factors (N), wave angles (v[/), and phase speed (c r ) vs. spanwise wavelengths (1) 
for a disturbance frequency of 2,000F1Z. (a) x/c = 0.25, (b) x/c = 0.45. 
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